function bestlambda=logitlev1(lambda,x,y) 

n=length(x);
r=length(lambda);
pred=cell(1,r);
Err=zeros(1,r);
for q=1:r
    l=lambda(q);
    pred{q}=zeros(n,3);
    for g=1:n
        sumpay=sum(x{g},2);
        pred{q}(g,:)=exp(l*sumpay);
        pred{q}(g,:)=pred{q}(g,:)./sum(pred{q}(g,:));
    end
    err=zeros(n,1);
    for i=1:n
        err(i)=-log(pred{q}(i,str2double(y(i))));
    end
    Err(q)=mean(err);
end

bestlambda=find(Err==min(Err));